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Abstract. In the present paper we study the equation 
of state of strange matter build of the u,d,s quarks in 
the framework of the MIT bag model. The scaling rela- 
tions with bag constant are discussed and applied to the 
determination of maximum frequency of a particle in sta- 
ble circular orbit around strange star. The highest QPO 
frequency of 1.33 kHz observed so far is consistent with 
the strange stars models for which the maximum QPO 
frequency is 1.7-2.4 kHz depending on the strange quark 
mass and the QCD coupling constant. The linear approx- 
imation of the equation of state is found and the parame- 
ters of this EOS are determined as a functions of strange 
quark mass, QCD coupling constant and bag constant. 
This approximation reproduces exact results within an er- 
ror of the order of 1% and can be used for the complete 
study of the properties of strange stars including micro- 
scopic stability of strange matter and determination of the 
total baryon number of the star. 
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1. Introduction 

The idea of the compact stars build of strange matter 
was presented by Witten ( 1984 ) and the models of stars 
were calculated using various models of strange matter by 
Haensel et al. ( |1986D and Alcock et al. ( [19861) . The main 
idea is that the u,d,s matter is ground state of matter at 
zero pressure (self-bound strange quark matter) i.e.: 



fio = ^i{P = 0) < M(^''Fe) = 930.4 MeV 



(1) 



There were two aspects of rather extensive studies 
of strange stars properties. The first one in the con- 
text of maximum rotational velocity of the star (Frieman 
& Olinto |1989[ , G lcnde nning [1989| , Prakash et al. |199C , 
Zdunik & Haensel 1990) was related to the announcement 



reporting the detection of the half millisecond pulsar in 
1989 (which after one year turned out to be erroneous). 
This observation stimulated detailed investigations of the 
limits on the rotational frequency of dense stars which 



excluded nearly all neutron stars models leaving strange 
stars as a possible explanation. 

At present the increasing interest in strange stars is 
a result of QPOs observations in low-mass X-ray binaries 
LMBX and some difficulties with the explanations of its 
properties by neutron star models under the assumption 
that QPOs represent Keplerian frequencies of the particles 
in an accretion disk (Kaaret et al. 1997, Kluzniak 1998, 
Zhang et al. 1998, Miller et al. 1998, Thampan et al. 1999, 
Schaab & Weigel 1999). In this paper the maximum val- 
ues of these frequencies are found for a broad set of the 
parameters describing a strange matter EOS. Strange star 
models are consistent with the maximum observed QPO 
frequency 1.33 kHz in 4U 0614-fG9 (van Straaten et al. 

pooo| ). 



Frieman & Olinto (1989) mentioned the approximation 



of the EOS of strange matter by the linear function (see 
also Prakash et al. 199C) . I present here the linear form 



of the EOS of strange matter with parameters expressed 
as a polynomial functions of the physical constants de- 
scribing MIT bag model (strange quark mass and QCD 
coupling constant). These approximate formulae allow us 
to write down the pressure vs. density dependence in a 
simple form P = a ■ {p — po). The consequence of this 
form are scaling laws of all stellar quantities with appro- 
priate powers of po- This linear EOS is complete in the 
sense that it contains not only pressure and energy den- 
sity, which enter the relativistic stress-energy tensor and 
are sufficient for the determination of the main parame- 
ters of the star (mass, radius), but also the formula for 
baryon number density (or equivalently baryon chemical 
potential), necessary in the microscopic stability consid- 
erations and determination of the total baryon number of 
the star. 

2. Strange stars and mELximum QPO frequency 

We describe the quark matter using the phenomenological 



MIT bag model (see, e.g., Baym, 1978). The quark matter 



is the mixture of the massless u and d quarks, electrons 
and massive s quarks. The model is described in detail 
in Far hi & Jaffc ( |l984|) , where the formulae for physical 
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parameters of a strange matter are also presented. There 
are the following physical quantities entering this model: 
B - the bag constant, ac - the QCD coupling constant 
and THs - the mass of the strange quark. It is necessary 
to introduce also the parameter pN - the renormalization 



point. Following Farhi & Jaffe (1984) we choose pN = 
313 MeV. 

The consequence of this model of strange matter is 
scaling of all thermodynamic functions with some powers 
of B. Knowing the EOS for given ac, rris and Bq we can 
obtain thermodynamic quantities for other value B from 
the following formulae: 



Pr 



IB] 



P[B] - P[Bo\ ■ (B/Bq), 



Bo 
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where the resulting EOS for B is determined for the same 
value of ac but for the strange quark mass given by the 
relation: 

ms{B)=m,{Bo)-{B/BoY/^ (3) 

The advantage of this approximations is the elimina- 
tion of one parameter {B) from the calculations of the 
EOS. The dependence of all parameters on B is very well 
defined and one can take this into account using simple 
scaling laws. It is therefore sufficient to study the EOS in 
two parameter space: ac and rris for chosen value of B^ 
and then obtain results for other B from Eqs (H). 

Strictly speaking the scaling relations (||) hold exactly 
for the parameter p^r rescaled by the relation Pn{B) = 
Pn[Bq) ■ {B/ BqY/^. However the renormalization point is 
fixed and this is the source of small discrepancy between 
the true results for different values of bag constant and the 
scaling relations presented in this paper. This difference 
is presented in Fig. |l] in the case of the baryon chemical 
potential at zero pressure calculated directly and using 
scaling relation ^ . We choose p,Q because this thermody- 
namic function is crucial in microscopic stability consid- 
erations. The models of strange quark matter correspond 
to two values of bag constant: JSq = 60 MeVfm"^ and 
i? = 90 MeVfm"^. The strange quark masses in this two 
cases fulfill the relation (||), namely TOs(60) = 200 MeV 
and TOs(90) = 221.33 MeV. 

We see that scaling formulae give exact results for ac = 
because then p^ does not enter the EOS. The maximum 
error of the scaling formulae (||) is less than 0.5%. 

From the microscopic point of view the strange quark 
matter is unstable at zero pressure for some values of ac 
and TOs discussed in this paper (large ac and nis), i.e. 
equation of state does not fulfill the relation (P . But from 
presented results we can obtain, using scaling formulae, 
the configurations which correspond to smaller value of 
B and are microscopically stable. The maximum value of 
B for which the strange quark matter is stable can be 
obtained from the equation: 

/ o \ 1/4 

po[Bo,ac,msiB^^y,/Bo)-^/^] [^^] = 930.4MeV(4) 



Fig. 1. The comparison between exact values of po, M 
and R obtained by the calculation of the EOS for B = 
60 MeV fm~'^ and B = 90 MeV fm"-^ and scaling relations. 
The values for B — 90 was rescaled to the value of Bq — 
60 using the relations /i — > /i • (B/i^o)"^^'*, rng —> ■ 
(S/Bo)-l/^ M M-[B/BoY'^ and R ^ R-{B/Bo)^'^. 
In the case of exact scaling relations the results for S = 60 
and rescaled results for B — 90 would be the same. 

The strange star configurations are calculated by solv- 
ing Oppenheimer-Volkoff equations in the case of spheri- 
cal symmetry. For given central pressure Pccntr we obtain 
stellar parameters such as gravitational mass A/, baryon 
number A, moment of inertia for slow rigid rotation / and 
gravitational surface redshift z. All these quantities are 
subject to the scaling formulae similar to those describing 
the EOS (Eq. ||). If we calculate the star for B^ with cen- 
tral pressure equal to Pccnti [5o] we know all parameters 
of the corresponding star with Pcontr = ^contr[-Bo] B / Bq in 
the model with bag constant B (cf. Witten 1984, Haensel 
et al. 1986). In general for the stellar parameter X the 



following equality holds: 



X[B,ac,ms] - X[Bo,ac,ms{B/Bo)-'^^] (B/Bo) 



(5) 



where /c = 1/2 in the case of mass and radius {X = R, M), 
k = 3/4 for X = A,k = 3/2 for / and fc = for z. 

In Fig. 1^ we see that scaling formulae (|^) reproduce 
the stellar parameters M vs. R with an error less than 
0.3%. 

These scaling laws so far have been mainly exploited 
in the case of maximum mass of the star (or equivalently 
maximum rotational frequency) . Of course these relations 
refer to all stellar configurations (e.g. the curve M{R)) 
and as an interesting example we can consider the point 
defined by the equality: 



GM 



R — Rms — 3Rg —6-2 



(6) 



which corresponds to the maximum possible frequency of 
a particle in stable circular orbit z/max around nonrotating 
star. Because the scaling properties of R and M are the 
same {'^ B^^^^) one can solve the Eq. ^independently of 
B. For stellar masses higher than the solution of the Eq. || 
the maximum Keplerian frequency of an orbiting particle 
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is defined by the marginally stable orbit at R^g and for 
lower by the stellar radius R (Fig. |^a) . In both cases this 
frequency scales as S^/^ and for other values of B the 
patterns of Fig. ^ do not change, provided one rescales 
the axes and uis (Eqs ^, ||, |^) . 

The maximum values of as a function of strange 
matter parameters are presented in Fig. ^3 for Bq — 
60 MeVfm^'^. For other values of B these results scale 
as ^1/2 i.e.: 
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which allows us to determine the absolute maximum value 
of i^max consistent with the requirement of the stability of 
strange matter at zero pressure (Eq. ^) by putting the 
maximum allowed value of B into Eq. (^. The result of 
this procedure vmax = i^max(.B,nax) is presented in Fig. 
We see that this limit is well above the highest observed 
QPO frequency i^obs — 1-33 kHz recently reported by van 
Straaten et al. ( 2000 ). Thus at present the observations 
of the highest QPO frequency do not constrain strange 
matter models unless we do not make the assumptions 
about the mass of the star in LMBX. To set some limits 
on strange matter parameters the QPOs must be observed 
at frequencies larger than 1.8 kHz. Our conclusion is true 
also for moderate rotation rates of strange star due to 
the initial increase of the frequency at marginally stable 
orbit which in the case of stars with mass slightly above 
the solution of Eq. ^ results in j/max very close to the value 
corresponding to the nonrotating stars (Zdunik et al. 2000, 
Datta et al. |200C| ). 

It should be mentioned that in principle analysis of the 
observational data can strongly support the existence of 
the marginally stable orbit. Such a c onclu sion have been 
recently presented by Kaaret et al. ( 199£ ) in the case of 
4U 1820-30. This interpretation is equivalent to the condi- 
tion R < Rms setting the lower bound for the mass of the 
star. This minimum mass limits are presented in Fig. ^ 
as a function of m^. and ac- Below this mass the innermost 
stable orbit is located at the surface of the star. 



3. Linear interpolation of EOS 

In the present section we discuss the interpolation of equa- 
tion of state of strange quark matter by the linear function 
of p. For a broad set of parameters ac and the depen- 
dence P{p) can be very well approximated by the linear 
function, being exactly linear at the rrig = limit. The 
EOS for strange quark matter is presented in Fig. ^. 

The main parameters of the matter are described by 
the following formulae: 



P{p) 
n{P) 



^(1 



no ■ 



efit)(p - 
1+14- 



3£fit. 

1 + £fit y poc' 



p 



-, 3/(4+efit) 



(8) 
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Fig. 2. The value of the frequency of a particle in the 
innermost stable circular orbit around nonrotating strange 
star for strange matter models with to^ = and Ws = 
200 MeV (panel a) and the maxima of those functions 
as a function of rus for ac — 0, 0.2, 0.4, 0.6 (panel b). 
The value of bag constant is fixed and equal to Bq = 
60 MeV fm-3 
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Fig. 3. The absolute maximum value of the frequency of a 
particle in stable circular orbit around nonrotating strange 
star vmax obtained for B = Bmax(TOs,ac) (Eq. ^ (panel 
a) and the corresponding values of stellar mass M - the 
minimum mass for which the existence of the marginally 
stable orbit is possible (panel b). 



where po and no are energy and number density of the 
strange quark matter at zero pressure. The second equa- 
tion is implied by the first law of thermodynamics. 

The equations (^) allow us to determine all thermody- 
namic quantities characterizing strange quark matter in 
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p[10"' g cm-=] 

Fig. 4. The linear approximation of the EOS for strange 
quark matter for three choices of {ms,ac): (0,0), 
(100,0.4), (250,0.6) (from left to right). The points rep- 
resent exact results obtained by the numerical calculation 
of the EOS for B = 60 MeVfm"^ and lines correspond to 
the approximate formula 



linear interpolation of the EOS, e.g. the baryon chemical 
potential: 



4 + £fit P 
1 + £fit Poc^ 



(l+efit)/(4+£fit) 



(9) 



where /zq — pqc^ /uq. 

Because for all thermodynamic quantities the scaling 
relations with B hold we can calculate all parameters 
needed for the linear EOS for chosen value of B and then 
rescale them using equations Thus the main point 
is to determine three parameters entering equations (^, 
as a functions of etc and mg: po^acrris), no(ac,r7is) and 
£fit (q^c J "T's)- These functions can be very well approxi- 
mated by the polynomial of with coefficients depending 
on the value of ac- The formulae obtained by the least- 
squares fit to the exact results are following: 



no 



)B 



■3/4 



Mo 



(oq + "^Loo + a's ™sioo 
= 0.28660Cy'^ 

(0.010788 + 0.0032046ac)/^^^^ 
-0.0044248a. 

« + a2"^siQ0 + a3"^sioo)^ 



(10) 



(11) 



1/4 



837.260/Q 
(46.616- 16. 848/Ca)/B 



1/2 
j3/4 



-fit 



(-10.482 -f4.5211/Ca)/B 

i'^slOO + £^3'^sl00 

(-0.035745 + 0.013366ac -t- 0.023246a^)/:B 
(0.0055488 - 0.0054232ac - 0.0069193a2)/;B 



(12) 

1/2 
■3/4 



where TOsIoo is strange quark mass in units 100 MeV 
i.e. 771,100 = m,c^ [MeV]/100, B = B [McVfm-^J/eo and 
Ca — I — 2ac/7r. 
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Fig. 5. The comparison between exact values of the main 
parameters of the strange matter (baryon chemical po- 
tential, particle number density, mass-energy density) at 
zero pressure and the parameter egt entering the linear 
approximate EOS of strange quark matter P = 1/3(1 -I- 
£fit)(p — Po) obtained by the interpolation of the EOS for 
B = 60 MeVfm^'^ (points) and approximate formulae ( [l0| - 
|l|) (fines). 

The energy-matter density at zero pressure is equal to 
Po = no^o/c^. The functions po{ms) and £fit(77is) entering 
directly Eq. ^ and p,Q{ms) for chosen values of ac and ac- 
curacy of the above interpolations are presented in Fig. ^ 

One can test the accuracy of the interpolation of the 
EOS by the linear function calculating the stellar con- 
figurations build of such a matter and comparing results 
with stars calculated using exact equation of state. Such 
a comparison is presented in Fig. ^. The relative error in 
this procedure is always less than 1% in the case of mass, 
radius and moment of inertia. 

Except for the case of strange matter with massless 
quarks, studied in detail in the literature, the linear form 
of the equation of state of self-bound matter have been 
used so far mainly for the stellar mass and radius determi- 



nation (e.g. Haensel & Zdunik 1989, Lattimer et al. 1990). 
The authors exploited the relation P = a{p — po), because 
the quantities sufficient in these considerations are pres- 
sure and energy density, which enter the stress-energy ten- 
sor and TOV equations of hydrostatic equilibrium. How- 
ever for the complete study of the properties of star the full 
microscopic description of the matter is needed, including 
baryon chemical potential and particle number density. 
In our EOS the formula for 77(P) (Eq. ^ enables us to 
determine the total baryon number of the stellar configu- 
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Fig. 6. The mass-radius relations for strange stars cal- 
culated using exact EOS (points) and the linear approx- 
imation of the EOS P ^ [p — pq) (lines) for B = 
60 MeVfm^'^. The curves correspond to four different 
parameters of strange matter (tos,Q!c): (0,0), (150,0.3), 
(200,0), (250,0.6) (from right to left). 

ration and use this model of strange matter, for example, 
for the discussion of the conversion of neutron stars into 



strang e star s (Olinto 1987 , Cheng & Dai 1996 , Bombaci & 
Datta |2000| ). During this process the total baryon number 
of the star is fixed. Using the relation (O) one can find 



the ranges of the values of rus, ac and B consistent with 
the requirement that strange matter is ground state of 
matter at zero pressure (Eq. |l|). As a result the allowable 
parameters of strange stars (e.g. M, vmax) can be found 
without complicated and time-consuming calculations of 
the microscopic properties of strange matter. 



4. Discussion and conclusions 

In the present paper I analyzed the accuracy of the scaling 
properties of strange matter parameters with the value of 
the bag constant B in the framework of MIT bag model 
of quark matter. The scaling formulae reproduce the main 
stellar parameters with an error less than 1%. This allow 
us to use them to determine the maximum possible fre- 
quency of a particle in stable circular orbit around strange 
star. The absolute maximum QPO frequency that can be 
accommodated by strange stars ranges from 1.7 to 2.4 kHz 
depending on the values of a strange matter parameters: 
strange quark mass and QCD coupling constant. Thus the 
present status of observational data — 1.33 kHz) 

cannot exclude strange stars as source of QPOs in LMBX. 
The frequencies of QPO in the very high range (larger than 
1.8 kHz), if observed, may set some bounds on parame- 
ters of strange matter, excluding large rUs and large 
We should mention that such a high values of vqpo can- 
not be understood in terms of neutron stars (Thampan et 
al. 



P = a{p — po) for which such a scaling properties are well 
known and exact. For this linear EOS the parameters of 
the strange stars scale with the powers of po for fixed value 
of a (analogous to the scaling laws with B, Eqs. 0). It 
should be stressed that these scaling properties are valid 
not only for static stellar configurations but also in some 
dynamical problems (e.g. parameters of the rotating star, 
Gourgoulhon et al. 1999). For fixed po one can apply also 



1999) 



The high accuracy of the simple scaling laws with B 
is strictly connected with the possibility of the approx- 
imation of the equation of state by the linear function 



the approximate scalings of the stellar parameters at the 
maximum mass p oint with appropriate powers of a (Lat- 
timer et al., |l990|) which have been recently confirmed by 
Stergioulas et al. (|l99E ) for stars rotating at Keplerian 
frequency. 

The presented linear approximation of the equation 
of state allow us to determine the dependence of many 
properties of strange stars on the physical parameters of 
a matter (tos, Q^c) using very simple form of the EOS. For 
a broad set of rUs and ac the values of Sat and po entering 
linear EOS can be very accurately obtained by a poly- 
nomial formulae. In particular the formula for the baryon 
chemical potential enables us to study the microscopic sta- 
bility of strange matter and make a complete discussion 
of the resulting constraints on strange stars parameters. 
The expression for baryon number density is necessary in 
the consideration of the conversions of neutron stars into 
strange stars with total baryon number conserved. 

It is worth noticing that the linear approximation of 
P{p) dependence (Eq. ||) can be used also for other mod- 
els of strange matter, self-bound at high density po. The 
expressions for n{P) and p{P) allow us to determine all 
microscopic properties of the matter given the values of 
Po and uq. 
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